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Abstract. We extend Halphen's theorem which characterizes the solutions 
of certain nth-order differential equations with rational coefficients and mero- 
morphic fundamental systems to a first-order nxn system of differential equa- 
tions. As an application of this circle of ideas we consider stationary rational 
algcbro-geometric solutions of the KdV hierarchy and illustrate some of the 
connections with completely integrable models of the Calogero-Moser-type. In 
particular, our treatment recovers the complete characterization of the isospec- 
tral class of such rational KdV solutions in terms of a precise description of 
the Airault-McKean-Moser locus of their poles. 



1. Introduction 

The purpose of this paper is twofold. First we prove an extension of Halphen's 
theorem, which characterizes the fundamental system of solutions of certain nth- 
order ordinary differential equations with rational coefficients to first-order nxn 
systems. In the second part of this paper we show how to apply Halphen's theorem 
to completely integrable systems of the Calogero-Moser-type, recovering a complete 
characterization of the isospectral class of all algebro-geomctric rational solutions 
of the KdV hierarchy. 

We start by describing Halphen's original result. Consider the following nth-order 
differential equation 

q n (z)y< n \z) + q n ^ 1 (z)y^ 1 Hz) + ■■■ + q (z)y(z) = 0, (1.1) 

where Qj{z) are polynomials, and the order of q n (z) is at least the order of Qj{z) 
for all Q < j < (n - 1), that is, 

q m (z) are polynomials, < m < n, (1.2a) 
q m (z)/q n {z) are bounded near oo for all < m < n — 1. (l-2b) 

Then the zeros of q n (z) are the possible singularities of solutions of (1.1). 

Assuming the fundamental system of solutions of (1.1) to be meromorphic, the 
following theorem due to Halphen holds. 
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Theorem 1.1. (Halphen [21], Ince [22, p. 372-375]) Assume (1.2) and suppose 
(1.1) has a meromorphic fundamental system of solutions. Then the general solution 
of (1.1) is of the form 

n 

y(z) = c m r m (z)e x ™ z , (1.3) 

m— 1 

where r m (z) are rational functions of z, \. m G C, 1 < m < n, and c m , 1 < m < n 
are arbitrary complex constants. 

Moreover, the converse of Halphen's theorem holds as well. 

Theorem 1.2. (Ince [22, p. 374-375]) Suppose r. m (z) are rational functions of z 
and A m ,c m G C, 1 < m < n. If ri(z)e XlZ , . . . ,r„(z)e A " z are linearly independent, 
then 

n 

y(z) = Y c m r m {z)e XmZ (1.4) 

m— 1 

is the general solution of an nth-order equation of the type (1.1), whose coefficients 
satisfy (1.2). 

Remark 1.3. We note that Halphen's main idea of proof in [21] consists of replac- 
ing the rational coefficients in (1.1) by appropriate elliptic coefficients (as discussed 
in [20]) followed by an application of Picard's theorem (cf., e.g., [22, p. 375-378]). A 
closer examination of his argument seems to reveal a lack of proof of the crucial fact 
that the associated differential equation with elliptic coefficients necessarily has a 
meromorphic fundamental system of solutions. A proof of Theorem 1.1 (and The- 
orem 1.2), using a completely different strategy, is provided in Ince's monograph 
[22, p. 372-375]. 

One of the principal aims of this note is to prove a first-order n x n system gener- 
alization of Halphen's Theorem 1.1 and its converse, Theorem 1.2, in Section 2. 

Analogous results hold for nth-order equations and first-order systems with periodic 
and elliptic coefficients. For a glimpse at the vast literature in these cases and their 
applications to completely integrable systems we refer the interested reader to [15]- 
[18], [38], [39] and the literature therein. 

In Section 3 we then apply Halphen's theorem to the problem of characterizing the 
isospectral class of all stationary rational KdV solutions. All such (nonconstant) 
solutions q are well-known to be necessarily of the form 

M 

q(z)= qoo -J2s e (s e + l)(z-( e )- 2 (1.5) 
t=\ 

for some q x G C, {Q}i<i<m C C, (J ^ Q for £' ^ £, and 

M 

s£ G N, 1 < I < M with ^ s e {s e + 1) = g(g + 1) (1.6) 

t=\ 

for some g G N, and the underlying spectral curve is then of the especially simple 
rational type 

y 2 = (E- qoo f a+1 . (1.7) 
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On the other hand, not every q of the type (1.5), (1.6) is an algebro-geometric 
solution of the KdV hierarchy. In general, the points Q must satisfy a set of 
intricate constraints. In fact, necessary and sufficient conditions on Q for q in (1.5) 
to be a rational KdV solution are given by 

Ef?TW = for fc = l, and*=l,...,M. (1.8) 

This result was first derived by Duistermaat and Griinbaum [10] (cf. p. 199) in 1986, 
as a by-product of their investigations of bispectral pairs of differential operators. 
We will provide an elementary derivation of this result on the basis of Halphen's 
theorem and an explicit Frobenius-type analysis in Section 3. 

For a fixed g G N, (1.6) and (1.8) yield a complete parametrization of all rational 
KdV solutions belonging to the spectral curve (1.7). In other words, they provide 
a complete characterization of the isospectral class of KdV solutions corresponding 
to (1.7). The constraints (1.8) represent the proper generalization of the locus of 
poles introduced by Airault, McKean, and Moser [5] in the sense that they explicitly 
describe the situation where poles are permitted to collide (i.e., where some of the 
s e > 1). 



2. Halphen's theorem for first-order systems 

This section is devoted to a generalization of Halphen's theorem (and its converse) 
to first-order systems. 

We briefly describe some of the notation used in this section. /„ denotes the identity 
in C™. An to x to diagonal matrix D = {dj6j > k)i<j,k<m wu l occasionally be denoted 
by diag(c?i, . . . , d m ). The operation of transposition is denoted by the superscript 
t. Moreover, it will be convenient to denote the set of all m x n matrices whose 
entries are rational functions with respect to z G C by TZ mxn , the subset of TZ mxn 
with rational entries bounded at infinity by 7\L™ X ™. 

We recall that for T G 7^™ x " invertible and differentiable with respect to z, the 
transformation y(z) = T(z)u(z) turns the first-order system of differential equations 
y'(z) = A(z)y(z) into the system u'(z) = B(z)u(z), with B{z) = T(z)- l (A{z)T{z)- 
T(z)). 

Definition 2.1. (i) Two matrices A,B G JZ nxn are called of the same kind if there 
exists an invertible matrix T G TZ nxn such that 

B(z) = T(z)-HA(z)T(z)-T'{z)). (2.1) 

(ii) B G TZ nxn is called reduced of order k if Bj^ = Sj+ij for all 1 < j < k and 
l<£<n. 

Our approach, including the notion of matrices being "of the same kind", was 
inspired by Loewy [29]. The relation of being of the same kind is obviously an 
equivalence relation on lZ n x ™ . The relation of being of the same kind is obviously 
an equivalence relation on lZ nxn . 
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Lemma 2.2. Suppose that A G TV^ n is reduced of order k — 1. Then either 
Ak,k+i — ■ ■ ■ — Ak, n — 0, or else there exists a matrix B G W£ n of the same 
kind as A and also reduced of order k — 1 but with the additional property that 
Bk,k+i(oo) ^ 0. Moreover, A(oo) and B(oo) have the same eigenvalues counting 
algebraic multiplicities. 



Proof. We assume not all of the entries Ak,k+i,- ■ ■ , Ak t n are equal to zero. Consider 
the (n — k + 1) x (n — k) matrix in the lower right corner of A and denote it by 
R. Suppose that r is the largest nonnegative integer such that z r R\^{z) remains 
bounded near infinity for every j G {1, ...,n— k}. Then there exists an £ £ 
{1, . . . , n — k} such that z r R\^(z) does not vanish at infinity. Denote the constant 
(n — k) x (n — k) matrix, which achieves the exchange of columns 1 and I of R(z), 
by C. Then the first row of z r R(z)C is bounded at infinity and the first entry in 
that row does not vanish at infinity. Next, define 

r(*)=fn JU, (2-2) 



x z r C / 

where Ik is the k x k identity matrix. Let 

- (few fcS) ■ < 2 - 3 ' 

where j4i ; i(^) and Ai^,{z) are square matrices with k and n — k rows, respectively. 
Then 

Since only the last row of ^1.2(2) is different from zero, and since that row equals 
the first row of R(z), the matrix T(z)~ 1 A(z)T(z) remains bounded at infinity and 
its first k — 1 rows are the same as those of A(z). The matrix C was chosen so that 
the first entry in the last row of z r Ai t2 {z)C does not vanish at infinity. Since 

lim T(z)" 1 T'(z) = 0, (2.5) 

z — »oc 

we conclude that B = T- 1 (AT - T') G W£ n is reduced of order k - 1 and that 
B k ,k+i{oo) ^ 0. 

Finally we prove that A(oo) and B(oo) have the same eigenvalues counting algebraic 
multiplicities. Since T^oo) might not exist, we first compute 

dct( lim ((rMT)(z) - A/J) = lim dct((T- 1 AT)(z) - A/„) 

Z — >OQ Z—*OD 

= lim det(A(z) - \I n ) = dct( lim (A(z) - XI„)). (2.6) 

z — ^00 z — >oc 

By (2.5), the left-hand side of (2.6) is the characteristic polynomial of B(oo), while 
the right-hand side is the characteristic polynomial of A(oo). This completes the 
proof. □ 

Lemma 2.3. Assume that A G TZ^ n is reduced of order k — 1 and suppose that 
Ak,k+i(oo) ^ 0. Then there exists a matrix B G TV^ n of the same kind as A, which 
is reduced of order k. Moreover, A(oo) and B(oo) are similar and hence isospectral 
(i.e., their eigenvalues, including algebraic and geometric multiplicities, coincide). 
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Proof. Let T £ 7^«x« denote the n x n matrix obtained from the identity matrix 
/„ by replacing its (k + l)st row by 

{—Ak,i, • • • ' ~Ak,k, 1, -A k ,k+2, ■ ■■ , -Ak, n )/Ak,k+i- (2-7) 

T _1 is then the matrix obtained from the identity matrix I n by replacing the 
(k + l)st row by (A^i, . . . , A k _ n ). Note that the entries of T and T~ x are ratio- 
nal and bounded at infinity. Hence the matrix B = T~ 1 (AT — T") has rational 
entries bounded at infinity. A straightforward calculation then shows that the 
first k rows of B have the desired form. Since T and T 1 are bounded at infinity, 
lim^oo T(z)- 1 T'(z) = and hence B(oo) = T(oo)- 1 A(oo)T(oo). □ 

Theorem 2.4. Let Q 6 TZ^ n and suppose that the first-order system y'(z) = 
Q(z)y(z) has a meromorphic fundamental system of solutions. Then y' — Qy has 
a fundamental matrix of the type 

Y(z) = R(z) cxp(diag(Ai2;, . . . , \ n z)), (2.8) 

where \\, . . . , A„ are the eigenvalues of Q(oo) and R € lZ nxn . 

Proof. The theorem will be proved by induction on n. Let n = 1. Any pole of Q(z) 
must be of first-order with an integer residue, that is, 

N 

Q(z) = A 1+ £^^, (2.9) 

with mi, . . . , mjy e Z. Then Y(z) = Ilfcit 2 — a i) me exp(Aiz) proves the claim for 
n= 1. 

Next, let n be any natural number and assume that Theorem 2.4 has been proven 
for any natural number strictly less than n. 

By hypothesis, Q E TZ^ n and Q(z) can be regarded to be reduced at least of 
order zero. We denote the eigenvalues of Q(oo) by Ai, . . . , A„. Repeated, perhaps 
alternating, applications of Lemmas 2.2 and 2.3 then yield the existence of an 
integer k G {1, . . . , n}, a k x k matrix B\(z), an (n — k) x k matrix B 3 (z), and an 
(n — k) x (n — k) matrix B±{z), such that 

has the following properties: 

(1) B e TZ n J n . 

(2) B is of the same kind as Q, that is, there exists an invertible matrix T 6 
U nxn guch that B ( z j = T{z)- 1 {Q{z)T{z) - T'{z)). 

(3) 5i(z) is reduced of order k — 1. 

(4) After a suitable relabeling of the eigenvalues of Q(oo) the eigenvalues of 
-Bi(oo) are Ai, . . . , Afc and the eigenvalues of -64(00) are Afc+i, . . . , A„. 

(5) The first-order system u'(z) = B(z)u(z) has a meromorphic fundamental 
system of solutions with respect to z G C. 

We now have to distinguish whether k = n or k < n. In the case k = n, B(z) = 
Bi(z), and the system u'(z) — B(z)u(z) is equivalent to the scalar equation 

u^\z) = B nil (z) Ul + ■■■ + B ntn (z)u^- 1} (z). (2.11) 
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In this case Halphen's theorem, Theorem 1.1, and the relations Uk(z) = (z) 
and y{z) = T(z)u{z) prove our claim. 

Next, assume that k < n. If w is any solution of w'(z) = Bi(z)w(z), choose a 
solution v of the nonhomogeneous system 

v'(z) - B 4 (z)v{z) = B 3 (z)w{z). (2.12) 

Then u = (it;,?;)' is a solution of u'(z) — B(z)u(z) and hence meromorphic. Thus 
every solution of w'(z) = Bi(z)w(z) and, choosing w(z) = in (2.12), also every 
solution of v'(z) — B 4 (z)v(z) is meromorphic. By the induction hypothesis, there 
exist matrices i?i e K kxk and R 4 e ft(n-fc)x("-*0 such that 

Ui(z) = J Ri(z)diag(cxp(Aiz), . . . ,exp(A fe z)) (2.13) 

is a fundamental matrix of w'(z) = Bx{z)w{z) and 

U 4 (z) = R 4 (z) diag(cxp(A fe+ iz), . . . , cxp(A„z)) (2.14) 

is a fundamental matrix of v'(z) = B 4 (z)v(z). 
Next define 

U 3 (z) = U 4 (z) j Z d£Ur\QB 3 (<;)Ui(Q. (2.15) 

Then each column of 

<"«> 

is a solution of u'(z) = B(z)u(z) and U is indeed a fundamental matrix of u'(z) = 
B{z)u(z) since det(C/(z)) = det{U\(z)) &et(U 4 (z)) ^ 0. It remains to show that 

U 3 (z) = R 3 (z) diag(exp(Aiz), . . . , exp(A fc z)) (2.17) 

for some matrix R 3 G 7^("-fc)xfe rp^ e entry j n row j anc [ column I of the matrix 
U 4 X (£)B 3 (£)\J\{£) equals Pj,t{Q exp[(A£ — \k+j)(], where pj^ is a rational function, 
that is, 

N M M r , 

Pi AO - E <h*.rC + EE tfzty (2 - 18) 

r=0 r=l s=l ^ rJ 

for appropriate choices of the parameters <Lj^ r , bj^^ iS , and pairwisc distinct z r . 
Next we recall that 



d(Ce xc = f{s, A, z) (mod (c A *C(z))), (2.19) 



where 



f(s,\z) 



(0 ifs>0 

(-sl'i 1 )! Ei ( Az ) ifA^0ands<-l 

ln(z) if A = and s = — 1 

if A = and s < -2 



(2.20) 



and that the exponential integral Ei(-) has a logarithmic branch point at zero. 
Therefore, if Xi ^ Afc+j, 

{U 4 1 U 3 ) j , l {z) = f dC Pj M> {Xl ~ Xk+i)C (2-21) 
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M 



e^-^+i)*S jtt (z) + ^2 



c 



: (A,-A fe+3 ), r Ei((\t - \ k+j )(z - Zr)) 



for appropriate rational functions Sjj and appropriate complex numbers Cj,^ r . 
However, since the entries of U 3 {z) and U^{z) must be meromorphic, all of the 
numbers Cj t t, r must necessarily vanish. If Xg = Xk+j a similar conclusion shows that 
no logarithmic terms appear so that in either case (U^ Us)j t e(z) G e^ Xe ~ Xk +^ z C(z). 
Hence we obtain 



U^(z)U 3 (z) = diag(e- A '=+ lZ , . . . , e - A " z )S(z) diag(e AlZ , . . . , e AfcZ ), (2.22) 
where S G n (n - k '> xk is the matrix with entries S 3 j. Thus, R 3 = R^S G TZ (n -^ xk . 



Remark 2.5. If Q(oo) = 0, the transformation z = ^ immediately shows that 
( = 0, is a regular singular point of our differential equation. This implies that 
the fundamental system at £ = is of the form (cf., e.g., [36, Sect. 23]) Y(Q = 
U(()( m f° r ICI < Co, where U(() is a holomorphic matrix for \(\ < ( . Hence 
our fundamental system is meromorphic on the whole Riemann sphere and must 
therefore be a purely rational matrix. 

Next we consider two examples. 

Example 2.6. The first-order 2x2 system 



This seems to suggest consideration of even more general systems of the type 
z~ q Y'{z) = A(z)Y(z), with q > 0, rather than the case q = only. But Theo- 
rem 2.4 can not hold in general for q > as shown by the following elementary 
counterexample. 

Example 2.7. The first-order 2x2 system 



has no solution in terms of elementary functions, although it clearly has a mero- 
morphic fundamental system. The particular case m = 1 represents the well-known 
Airy equation. 

Remark 2.8. In the case where all eigenvalues 1 < Xj < n of Q(oo) arc distinct, we 
now sketch an alternative proof of Theorem 2.4, based on Theorem 12.3 in Wasow's 
monograph [37]. Since Theorem 12.3 in [37] only applies to appropriate sectors of 
the complex plane with vertex at the origin, we argue as follows. First one can 
find a sufficiently small sector S3, which does not contain any separation rays. (We 
recall that a ray (i.e., a half line), where Re(XjZ — Xkz) — for some pair of distinct 
integers j, k, is called a separation ray.) Then one chooses two other sectors Si, S2 
with opening angles <f)j < it, j = 1, 2, such that Si U S2 U S3 = C\{0}. It is then 
possible to show that the transition matrix from sector Si to sector S2 equals the 



□ 




is solved by 
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identity matrix. Hence, the solution of the form Y(z) = R(z) exp(diag(Ai, . . . , A„)z) 
in sector Si is valid in sector £2 too and thus can be continued into S3 since by 
hypothesis, the sector S3 contains no separation rays. 

Finally, we turn to a converse of Theorem 2.4. 

Theorem 2.9. Suppose R e TZ nxn , det(R) ^0, and Ai, ... , X„ G C. Then 

Y(z) = R(z) exp(diag(Aiz, . . . , X n z j) (2.23) 

is a fundamental matrix of a first-order linear system of differential equations 
y'(z) — Q(z)y(z), where Q G K" xn and Q(z) is of the same kind as a ma- 
trix in TV^ n . In fact, Q(z) is of the same kind as the constant diagonal matrix 
diag(Ai, . . . , A„). 

Proof. Since 

Q(z) = R{z) diag(Ai, . . . , A„)i?(z)- 1 + R'(z)R{z)-\ (2.24) 
we choose T(z) — R(z)^ 1 and hence obtain T" = —R~ 1 R'R~ 1 and thus, 

Q = T- 1 (diag(A 1 , A n )T - T'). (2.25) 

Hence, Q(z) is of the same kind as the constant matrix diag(Ai, . . . , A„). □ 

3. Some applications to rational solutions of the stationary KdV 

HIERARCHY 

In this section we describe the connections between the preceding results and 
infinite-dimensional completely integrable Hamiltonian systems. For reasons of 
brevity we will only consider the simplest case of the KdV hierarchy, and in ac- 
cordance with Sections 1, 2, only study its stationary rational solutions bounded 
at infinity (cf. [1], [3]-[7], [19], [23], [24]-[28], [30], [31], [32], [33], [35], [40] and the 
literature cited therein). The principal results on the stationary KdV hierarchy as 
needed in this section are summarized in the appendix, and we freely use these 
results and the notation established there in what follows. 

The rational KdV solutions bounded at infinity are usually discussed in a time- 
dependent setting and the dynamics of their poles is in an intimate relationship 
with completely integrable systems of the Calogero-Moser-type. In our discussion 
below, the time-dependence will generally be suppressed and only occasionally be 
mentioned in connection with particular isospectral deformations of rational solu- 
tions of the KdV hierarchy. Our principal focus will be on stationary (isospectral) 
aspects of these rational KdV solutions and the implications of Halphen's theorem 
in this context. 

We start by quoting a number of known results on stationary rational KdV solutions 
bounded at infinity. 



Theorem 3.1. Let WeN and {zj}i<j< N C C. 

(i) (Airault, McKean, and Moser [5]) Any rational solution q of (some, and hence 
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infinitely many equations of) the KdV hierarchy, or equivalently, any rational 
algebro- geometric potential q, is necessarily of the form 

N 

q(z) = qoo -2j2(z-z j r 2 , (3.1) 
for some qoo G C and with N G N of the special type N = g(g + l)/2 for some 

<?eN. 

(ii) {Airault, McKean, and Moser [5] {see also [38])) If one allows for "collisions" 
between the zj, that is, if the set {zj}i<j<N clusters into groups of points, then the 
corresponding rational algebro- geometric potential q is necessarily of the form 

M 

q(z)=q 00 -Y J Sl{si + l)(z-Q)- 2 , (3.2) 
1=1 

where for some g G N. 

{zj}i<j<N = {Q}i<£<M C C, with d pairwise distinct, (3.3a) 
Si G N, 1 < £ < M, 

M 

s i( s i + !) = 2N f or some N e N of the type N = g{g + l)/2. (3.3b) 

(=1 

(Hi) The extreme case of all zj colliding into one point, say £i, that is, {zj}i<j<N = 
{Ci} C C yields an algebro- geometric KdV potential of the elementary form 

q(z) = q 00 -g(g + l)(z-( 1 )- 2 , ff GN (3.4) 

and no additional constraints on Q\ G C. 

(iv) In all cases (i)-(iii), if q is a rational KdV potential (i.e., if g G N and the 
points Zj (resp. Q) satisfy appropriate restrictions, cf. Theorem 3.5), the underlying 
rational hyperelliptic curve JC g is of the especially simple form 

IC g :y 2 = (E- qoo )^+\ (3.5) 

In particular, the potentials (3.1), (3.2), and (3.4) are all isospectral (assuming 
(3.1) and (3.2) are algebro- geometric KdV potentials, of course). 

(v) (Weikard [38]) q is a rational KdV potential if and only if ip" + (q — E)ip = 
has a meromorphic fundamental solutions (w.r.t. z) for all values of the spectral 
parameter E G C. 

(vi) If q is a rational KdV potential of the form (3.2), then y" + qy = Ey has 
linearly independent solutions of the Baker- Akhiezer-type 

ME, z) = {± E^)- g ( ft ( ± £ 1/2 - Mz))y ±EU2z , (3.6) 

E G C\{ 9oo }, z G C, 
with Hj(z) — i>j(z) 2 , 1 < j < g the zeros of F g (z,x) as defined in (A. 12). 

(To avoid annoying case distinctions we will in almost all circumstances exclude 
the trivial case N = g = in this section.) 
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Remark 3.2. (i) It must be emphasized that for JV > 1, not any potential q of 
the type (3.1) is an algebro-geometric KdV potential. In fact, for N > 1, there 
exist nontrivial constraints on the set {zj}i<j<N for (3.1) to represent an algebro- 
geometric KdV potential. For instance, if the zj in (3.1) are pairwise distinct, then 
Airault, MacKean, and Moser [5] proved that 

N 1 

E (z ._ z .,) 3 =° for. 7 = l,...,iV (3.7) 

are necessary conditions for q in (3.1) to be a stationary KdV potential. In the 
case of collisions (i.e., if si a > 1 for some 1 < £q < M) the necessary constraints 
on {Q}i<£<m are more involved than in the nondegenerate case above and a com- 
plete description of all constraints were originally obtained by Duistcrmaat and 
Griinbaum [10] in 1986. An alternative proof of their result will be given in Theo- 
rem 3.5 below. 

(ii) In connection with Theorem 3.1 (ii) one might naively expect that any decom- 
position of g(g + 1) = X^=i s e( s £ + 1) can actually be realized for some choice 
°f {Q}i<t<M with Q ^ Q> for I ^ £'. However, the simple counterexample 
q(z] ) = -6(z-Cir 2 -6(z-C 2 r 2 , which satisfies KdV 3 (q) = -5670(Ci - C2) 2 (Ci + 
Q 2 — 2z)(z — (i)~ e (z — C2)~ 6 , quickly destroys such hopes. 

(iii) Strictly speaking, the version of Theorem 3.1 (v) proven in [38] assumes in 
addition to q being rational, that q is bounded at infinity. However, assuming that 

q(z) = az k + 0(z fe_1 ) for some a ^ and k G N, 
a simple inductive argument using (A.l) proves 




+ o(^' fe - 2 ), i>i, 



using the usual convention (for j = 1) that products over empty sets are put equal 
to one. Thus, since /j cannot vanish in this CMSC. el rational q unbounded at infinity 
cannot satisfy any of the stationary KdV equations (cf. (A. 8)). 

Before we discuss additional facts, we briefly pause and mention some of the ingre- 
dients entering the proof of items (i)-(v) in Theorem 3.1. We start with a fairly 
complete treatment of item (iii) and for simplicity of notation put goo = Ci — and 

q g (z) = -g(g+l)z- 2 , g e N, z e C\{0}. (3.8) 

From [2, Ch. 10] one infers that (E e C\{0}, z £ C) 

im*.*) = [t^^^^y^ (3.9) 

are linearly independent solutions of ip" + (q g — E)i[) = 0, E £ C\{0}. Thus, one 
concludes that 

g 

il>+(E,z)1>-(E,z) - n (l - ^) for some Kj £ C, 1 < j < g. (3.10) 
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Hence a comparison with (A.12)-(A.15), (A.19)-(A.24) yields 

9 

F g (E,z) = Y[(E-n j (z)), n j (z) = Kj z- 2 , l<j<g, (3.11) 

where F g (E,z) denotes the polynomial of degree g with respect to E associated 
with q g (z) in (3.8), as introduced in the appendix. Thus, q g (z) is a KdV poten- 
tial satisfying KdV g (q g ) = for a particular set of constants {ce}i<e< g in (A. 10). 
However, taking into account the simple form of q g (z) in (3.8), homogeneity con- 
siderations in connection with the corresponding fj and (A. 25) then yield in the 
special case q(z) = q g (z), 

ce = 0, 1 < t < g, (3.12) 

F g (E,z) = F g (E,z), f j (z) = f j (z), l<j<g, (3.13) 

fj(z) = djZ- 2j for some dj £ C\{0}, 1 < j < g, (3.14) 

f k+1 (z) = 0, s-KdV k (q g ) = 0, k>g, (3.15) 

y 2 = E 2g+1 , that is, E m = 0, < m < 2g (3.16) 

(and of course Co = fo{z) = fo(z) = 1). This yields item (iii) and part of item 

(iv) . Since q in (3.1) and (3.2) in the special case q^ = satisfies q(z) = 

\z\^oo 

2Nz~ 2 (l + 0(|z| -1 )), one infers that f k +i = for some fceN can only happen if 
N = k(k + l)/2 for some k e N. This illustrates N = g(g + l)/2 and (3.3b). Item 

(v) in [38] follows from a careful combination of Frobenius theory for second-order 
linear ordinary differential equations in the complex domain, Halphen's theorem, 
Theorem 1.1 (for n = 2), and some of the algebro-geometric formalism briefly 
sketched in the appendix. As a by-product of a proof of item (v) one shows that 
i>"(z) — cz~ 2 ip(z) — Etp(z), z s C\{0} has a meromorphic fundamental system of 
solutions for all E S C if and only if c 6 C is of the special form c = s(s + 1) for 
some s G No- This illustrates why collisions necessarily must happen as described 
in (3.3a). This fact was already known to Kruskal [28] in 1974. That q in (3.1), 
(3.2), and (3.4) are all isospectral KdV potentials, that is, they all belong to the 
same algebraic curve (3.5) (assuming (3.1) and (3.2) satisfy the additional restric- 
tions to make them algebro-geometric KdV potentials, of course) can be shown by 
several methods. Either by invoking time-dependent KdV flows as in [5], or by 
commutation techniques (i.e., Darboux-type transformations) as in [3], [11], [30], 
[31] (cf. also [14]). This fact also follows from the results in [38]. Finally, iden- 
tifying ip±(E, z)/ip±(E, zq) with the two branches of the Baker- Akhiezer function 
i>(P, z, zq), P = {E,y) in (A. 20), a combination of (A. 12), (A. 23), and the normal- 
izations 

lim tP ± (E 7 z)exp(TE 1/2 z) = 1, lim i/>±(E,z) exp(TE 1/2 z) = 1, 

|z|^oo |B|^oc 

then proves ip+(E, z)ip-(E, z) = ESF g (E,z) = Uj=i ^1-^^) , and hence (3.6). 

Finally, we study the precise restrictions on the set of poles {zj}\<j<N = {Ce}i<t<M 
for q in (3.2) to be a KdV potential. 
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Lemma 3.3. Suppose the function q has a Laurent expansion about the point zq G C 
of the type 

oo 

g (zH^(z-z r 2 , (3.17) 

where qo = —s(s + 1) and, without loss of generality, Rc(2s + 1) > 0. Define for 
a G C, 

fo(a) = -a(a - 1) - q Q = (s + a)(s + 1 - a), (3.18) 

<*{*)= Ufo(*+j), cj(a) = ^fft^M , jeN, (3.19) 

OO 

w(a,z)=J2cA<7)(z-zoy +j , (3.20) 

3=0 

v(a, z) = ^ {o -, z ) = jr / (^L + a log(z - z )) (z z y+i. (3.21) 

If 2s + 1 is not an integer, then y" + qy = /ias i/ie linearly independent solutions 
yi = w(s + 1, •) and y 2 = w{—s, ■). If 2s + 1 is an integer, then y" + qy = has the 
linearly independent solutions yi — w(s + 1, •) and y 2 = v(—s, •). 

Moreover, y" + qy = has a meromorphic fundamental system of solutions near z 
if and only if s e No and c 2s +i(— s) = 0. 

This is a classical result in ordinary differential equations (cf., e.g., [22], Chs. XV, 
XVI). A recent proof can be found in Section 3 of [38]. 

Definition 3.4. Let q be a rational function. Then q is called a Halphen potential 
if it is bounded near infinity and if y" + qy = Ey has a meromorphic fundamental 
system of solutions (w.r.t. z) for each value of the complex spectral parameter E G 
C. 



Of course every constant is a Halphen potential. Moreover, by Theorem 3.1 (v), q 
is a Halphen potential if and only if it is a rational KdV potential (i.e., if and only 
if it satisfies one and hence infinitely many of the equations of the stationary KdV 
hierarchy) . 

Theorem 3.5. Let q be a nonconstant rational function. Then q is a Halphen 
potential if and only if there are MeN, S(6N, 1 < £ < M, q^ e C, and pairwise 
distinct Q G C, I = 1, . . . , M, such that 

M 

q(z) =Qoo-J2 s ^ + - ^ ( 3 - 22 ) 
t=\ 

and 

E (T-CA^+i = ° for k = l,...,s e audi =1,...,M. (3.23) 

Moreover, q is a rational KdV potential if and only if q is of the type (3.22) and the 
constraints (3.23) hold. In particular, for fixed g, the constraints (3.23) characterize 
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the isospectral class of all rational KdV potentials associated with the curve y 2 — 
(E - qoo) 29+1 , where g(g + 1) = Y,?=i s ^ + 

Proof. By Theorem 3.1 (v), it suffices to prove the characterization of Halphen 
potentials. Suppose that q is a nonconstant Halphen potential. Then a pole zq of 
q is a regular singular point of y" + qy = Ey and hence 

oo 

q{z)-E^Y.^^-^Y- 2 
3=0 

in a sufficiently small neighborhood of zq, where Q 2 is a hrst order polynomial 
in E, while Qj for j ^ 2 are independent of E. The indices associated with zo, 
defined as the roots of 0(0 — 1) + Q = (hence they are ^-independent), must be 
distinct integers whose sum equals one. We denote them by — s and s + 1 where 
s > and note that Q = — s(s + 1). We intend to prove that Qij+\ = whenever 
j G {0, s} by applying Lemma 3.3. Proceeding by way of contradiction, we thus 
assume that for some nonnegative integer k G {0,...,s}, Q2ZC+1 7^ and k is the 
smallest such integer. 

We note that /o(- + j) are positive in (— s — 1, — s + 1) for j = 1, 2s, whereas 
/o(- + 2s + 1) has a simple zero at — s and its derivative is negative at — s. Next 
one defines 

2s+l 2s+l 

7o(<t)= I] hip + j) and 7i(<7)= I] M° + j)- ( 3 - 24 ) 

70 and 71 have simple zeros at — s and and 7 (— s) and 7^(— s) are negative. 

The functions Co = 70 and c\ — Q171 are polynomials with respect to E. Actually, 
Co has degree zero in E and c\ is constant but might equal zero. Hence the relations 
(3.25), (3.26), and (3.27) below are satisfied for j = 1. Next we assume that for 
some integer I <G {1, s}, the functions cq, c 2 e-i are polynomials in E and that 
the relations 

czj-iW) = 72 J - 2 ( ( r)Qr 1 + 0(&- 2 ), (3.25) 

r ( , /72 J -l( ( T)Q2 fe+ l0r^ 1 +O(^- 2 ), j-l>k, 

C2i - l(<7) = \0, j-Kk, (3 - 26) 

72 J -- 2 (-s)=72i-i(-s)=0, 72 J -2(-s),7 2j --i(-s)<0 (3.27) 

are satisfied for 1 < j < I. Using the recursion relation (3.19) we then obtain that 
c 2t {a) an d c 2£+i( cr ) are polynomials in E and that 

, , 121-2(0-) . n(T7t-i\ 
/o(<T + It) 



C2l+\{p) 



a )=\ /o(<r+2l+l) Q2k+lQ 2 +0(E ), £>k, 

0, £ < k. 



Letting = I21-2/ + 2£) and j 2l+1 = (72^-1 + 72(£-fc))//o(- + 2^+1) we find 
that the relations (3.25), (3.26), and (3.27) are satisfied for j = £ + 1. Hence an 
inductive argument proves that c 2s +i is a polynomial in E and that 

/ n 72s-l (oj + 72(s-fc)(c r ) ns _ fe n . ps _ fe _i, 

C2s+l((T)= / (<7 + 2 a + l) Q2fc+1 ° 2 + ° ( ^ } 
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= 72.s+iM02fc+ior fe + oiE*-*- 1 ). 

But both 72s-i+72(s-fc) an d /o(- + 2s + l) have simple zeros at — s so that 72 S +i(— s) 
is different from zero. Lemma 3.3 then shows that y" + qy = Ey has a solution 
which is not meromorphic whenever E is not a root of the polynomial C2 S +i(— s). 
This contradiction proves our assumption Q2fe+i i= wrong. 

Since Q\ — 0, we proved that if q is a Halphen potential with pairwise distinct 
poles Cij ■ • ■ 7 Cm, then the principal part of q about any pole Q is of the form 
—se(si + l)/(z — Q) 2 for an appropriate positive integer sg. Since q is bounded at 
infinity a partial fraction expansion then proves (3.22). This immediately implies 
that for zo = (e, 

This proves necessity of the conditions (3.22) and (3.23) for q to be a Halphen po- 
tential. To prove their sufficiency we now assume that (3.22) and (3.23) hold. Then, 
if zq denotes any of the points Q, one infers that the corresponding C2 Sf +i(— se) = 0. 
Lemma 3.3 then guarantees that all solutions of y" + qy = Ey are meromorphic 
and hence that q is a Halphen potential. □ 

Remark 3.6. (i) We emphasize again that the necessary and sufficient condi- 
tions on Q for q in (3.22) to be a rational KdV potential were first obtained by 
Duistermaat and Grunbaum [10] in their analysis of bispectral pairs of differential 
operators. Our approach based on Halphcn's theorem and a direct Frobenius-type 
analysis is a bit more streamlined since we aim directly at rational KdV solutions 
(and do not cover the case of the Airy equation) but there are undoubtedly some 
similarities in both approaches. 

(ii) We note that the restrictions (3.23) simplify in the absence of collisions, where 
se = 1, 1 < i < N. In this case (3.23) reduces to E^=i,j'^(^ - z j')~ 3 = °> 
1 < j < N, which represents the well-known locus introduced by Airault, McKean, 
and Moser [5]. This locus generated considerable interest, and especially its gener- 
alizations to elliptic KdV potentials and (elliptic) KP potentials, were intensively 
studied (cf., e.g., [4], [6], [7], [8], [12], [13], [23], [24]-[28], [32], [33], [34], [40]). The 
current derivation of (3.23) properly extends this locus to the case of collisions (i.e., 
to cases where some of the se > 1). Moreover, this appears to be the first system- 
atic derivation of this locus (with or without collisions) within a purely stationary 
approach (i.e., without involving special time-dependent KdV flows, etc.). 

(iii) For k = 1, conditions (3.23) coincide with the necessary conditions at colli- 
sion points found by Airault, McKean, and Moser [5] in their Remark 1 on p. 113. 
However, since there are additional necessary conditions in (3.23) corresponding to 
k > 2, this disproves the conjecture made at the end of the proof of their Remark 1. 

(iv) The genus g = 2 (N = 3) example, q 2 {z,t) = -6z(z 3 + 6t)(z 3 - 3i)~ 2 , t e C, 
with Zj = (Sty^ujj, ujj — exp(27rr//3), 1 < j < 3, explicitly illustrates the lo- 
cus in (3.23). One verifies that q 2 {t) satisfies the fcth stationary KdV equation, 
s-KdVfc(<?2(£)) = for all k > 2 and all t G R, as well as the 1st time-dependent 
KdV equation q 2lt = ^~ x q 2 , X xx + ^ 1 3q 2 q2,x (see, e.g., [4], [10]). 

Extensions of the stationary formalism described in this section to elliptic KdV 
potentials are in preparation. 
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The stationary KdV hierarchy 

In this section we review basic facts on the stationary KdV hierarchy. Since this 
material is well-known, we confine ourselves to a brief account. Assuming q to be 
meromorphic in C, consider the recursion relation 

/«,(*) = 1, f' j+1 (z) = 4-\/f (z) + qWf'jiz) + 2- 1 q'(z)f j (z) (A.l) 

for j G No (with / denoting differentiation with respect to z and No = N U {0}) and 
the associated differential expressions (Lax pair) 

dz 2 



L ? = ^2 ( A -2) 



3=0 

One can show that 



l?->, ge No. (A.3) 



P2g+1, L2 



= 2/: +1 (A.4) 



([■, •] the commutator symbol) and explicitly computes from (A.l), 

fo = 1, fx = 2" 1 g + ci, h = 8" V + S^iq 2 + + 02, etc., (A.5) 

where Cj G C are integration constants. Using the convention that the correspond- 
ing homogeneous quantities obtained by setting ce = for £ = 1, 2, . . . are denoted 
by fj, that is, 

fi=fi\c t =0,l<t<J> ( A - 6 ) 

one obtains 

3 

fj = J2cefj-e, 0<3<g- (A.7) 

The (homogeneous) stationary KdV hierarchy is then defined as the sequence of 
equations 

s-KdV 3 (<z) = 2f' g+1 = 0, g G N . (A.8) 

Explicitly, this yields 

s-KdVo(g) = q' = 0, s-KdVi(g) = 4"Y" + 2-^qq' = 0, etc. (A.9) 
The corresponding nonhomogeneous version of s-KdV s (g) = is then defined by 

£KdV„(«) = 2f' g+1 = 2 ]T Cg-j f'i+x = «. (A.10) 

where c = 1 and arbitrary complex constants. 

If one assigns to q^ = d e q/dz e the degree dcg(q^) = I + 2, / 6 No, then the 
homogeneous differential polynomial fj with respect to q turns out to have degree 
2j, that is, 

deg(/ j ) = 2i, j G No- (A.11) 
Next, introduce the polynomial F g {E 1 z) in E G C, 

F 9 (£,z) =Y,f g -i {z)& = Y[{E-h{z)). (A.12) 

3=0 3=1 
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Since f (z) = 1, 

- 2- 1 F>>(E, z)F g (E, z) + A-'F' g (E, zf + (E q(z))F g (E, zf 

= R 2g+ i(E,z) (A.13) 

is a monic polynomial in E of degree 2g + 1. However, equations (A.l) and (A. 10) 
imply that 

Z^F'g" -2{E-q)F g +q'F g =0 (A.14) 

and this shows that R2 g +i(E, z) is in fact independent of z. Hence it can be written 
as 

R2g+l{E) = ]J(E-E m ), {E m } < m <2g C C. (A.15) 
m=0 

By (A. 4) the nonhomogeneous KdV equation (A. 10) is equivalent to the commu- 
tativity of Li and P2g+i- This shows that 

[A s+ i,£ 2 ] = 0, (A.16) 

and therefore, if L>2tp = Eip, this implies that ffg+i^ = R-2g+i{E)ip. Thus 

[P2g+i,L2] = implies -P| 9 +i = R2g+i(L 2 ) by the Burchnall and Chaundy theo- 
rem. This illustrates the intimate connection between the stationary KdV equation 
f' g+1 = in (A. 10) and the compact (possibly singular) hyperelliptic curve JC g of 
(arithmetic) genus g obtained upon one-point compactification of the curve 

2g 

t g :y 2 = R2 9+1 {E)=Y[{E-E m ) (A.17) 

m— 

by joining the point at infinity, denoted by P^. Points P e fc g \{Poo} will be 
denoted by P — (E,y), moreover, the involution (hyperelliptic sheet exchange 
map) * on ft g is defined by 

*:£ g ^t g , P=(E,y)~P* = (E,-y),P£ 00 = P Too . (A.18) 

Introducing the meromorphic function </>(•, z) on JC g , 

<f>(P,z)=[y(P) + (l/2)F g (E,z)]/F g (E,z), P = (E,y) e t g (A.19) 
and the stationary Baker- Akhiezer function *p(-,z,z ) by 

V(P,^o)=expQ dz'<f>(P,z')j, PEJCgMPoo}, (A.20) 

one infers (for P = (E,y) e /C s \{Poo}, (z, z ) e C 2 ) 

L 2 iP(P,;z ) = Ei;(P,;z ), (A.21) 

P 2g+ ii>(P, ; zo) = yV(P, •, zo), (A.22) 

<1>(P, z, z )i>(P*,z, z ) - F g (E, z)/F g {E, z ), (A.23) 

Wty(P, ■, z ), 4>(P*, ; z )) = -2y(P)/F g (E, z ), (A.24) 

where W(f,g)(z) = f(z)g'(z) — f'(z)g(z) denotes the Wronskian of / and g. Thus, 
ip(P, z, Zo) and ^(P* , z, z a ) are linearly independent solutions of L 2 ip = Eip as 
long as E £ C\{E m }o< m <2 g . The two branches of tp{P, z, zo) will be denoted by 
tp±(E, z, zq), respectively. 
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The above formalism leads to the following standard definition. 

Definition .1. Any solution q of one of the stationary KdV equations (A. 10) is 
called an algebro-geometric KdV potential. 

For brevity of notation we will occasionally call such q simply KdV potentials. 
Finally, denoting E_ = (E , ■ ■ ■ , E 2g ), consider 
/ 29 / TP \ \ 1/2 00 

(nfi-T 1 )) =£-^ fe < 

V rn=0 V 7 7 fe=0 

1 N 

where c (E) = 1, a(E) = -- ^ E m , etc. 

Assuming that q satisfies the gth stationary (nonhomogeneous) KdV equation 
(A. 10), the integration constants ct in (A. 7) become a functional of the E m in 
the underlying curve (A. 17) and one verifies 

c e =c e (E), £ = 0,...,g. (A.25) 

Acknowledgment. We are indebted to Wolfgang Bulla for discussions on this 
subject. 
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